Coexistence of charge density waves, bond order waves and spin density waves in 

quasi-one dimensional charge transfer salts 
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Charge, spin, as well as lattice instabilities are investigated in isolated or weakly coupled chains of 
correlated electrons at quarter-filling. Our analysis is based on extended Hubbard models including 
nearest neighbor repulsion and Peierls coupling to lattice degrees of freedom. While treating the 
electronic quantum fluctuations exactly, the lattice structure is optimized self-consistently. We 
show that, generically, isolated chains undergo instabilities towards coexisting charge density waves 
(CDW) and bond order waves (BOW) insulating spin-gapped phases. The spin and charge gaps of 
the BOW-CDW phase are computed. In the presence of an interchain magnetic coupling spin density 
waves phases including a CDW or a BOW component are also found. Our results are discussed in 
the context of insulating charge transfer salts. 

PACS: 75.10.-b, 75.50.Ee, 71.27. +a, 75.40.Mg 



Quasi-one dimensional correlated electrons systems 
at quarter-filling show fascinating physical properties. 
A widely studied class of such materials are the so- 
called organic charge transfer salts like the Bechgaard 
salts (TMTSF) 2 X or their sulfur analogs (TMTTF) 2 X 
(X=PF 6 , AsF 6 ) These systems which consist of 

stacks of organics molecules forming weakly coupled one 
dimensional (ID) chains exhibit, at low temperature, a 
large variety of exotic phases such as triplet supercon- 
ducting [Q, spin density wave (SDW), charge density 
wave (CDW) and spin- Peierls (SP) phases f§. The com- 
pounds of the sulfur (TMTTF) 2 X series show strong 
charge localization at rather high temperature T p (as 
signalled e.g. by transport measurements) as well as a 
weak (to moderate) dimerization along the stacks Pl^). 
The insulating behavior observed in this regime has 
been interpreted as ID CDW fluctuations || or using 
a Hubbard chain Hamiltonian with an explicit dimeriza- 
tion [pi. At a significantly lower temperature T$p (typ- 
ically Tsp ~ 15K) a SP transition occurs together with 
a tetramerization along the chains Jl|,§|. In this system 
the interplay between charge ordering and lattice insta- 
bility is poorly understood. In particular whether the 
tetramerization (connected to the spin gap) is occurring 
simultaneously with a CDW transition is still controver- 
sial. So far the existing theory of the SP phase (7| do 
not consider the possibility of a coexisting 2kp CDW. 
Interestingly enough, CDW fluctuations were seen by X- 
ray diffuse scattering |f| in the SDW phase of the (more 
metallic) (TMTSF) 2 PF 6 compound. 

In this paper, our aim is to investigate by numeri- 
cal exact diagonalization (ED) techniques, the interplay 
between the electron repulsion and the electron-phonon 
coupling in the case of an adiabatic lattice. We focus 
on the competition or the cooperative behavior between 
charge ordering and lattice instabilities. Such a problem 
has been addressed in a number of previous studies Pfl0| 



where several interesting modulated phases have been 
proposed. However, so far, no systematic investigation 
of the full phase diagrams has been carried out. Indeed, 
the suggested translation symmetry broken states |l(]] 
were found in a restricted variational space. We shall 
here re-examine these issues in order to determine the 
absolute stability of the various competing phases. 

We use extended ID Hubbard models at quarter-filling 
(n = 1/2) coupled with some (classical) phonon field, 

Hid ='^2 t ( i ) {4;<T C i+U<r + h.c.) + U'^2 n »;T n *;l 

+ V niTll + 1 + H P h ' ( 1 ) 
i 

where n^ a — c[. a Ci- a and = rij ; j + n^j.. We have in- 
cluded a nearest neighbor (NN) interaction V as its role 
will become clear in the following. Local deformations 
of the molecules can produce changes of the on-site (or 
molecular) orbitals energies and can simply be taken into 
account by a Holstein term mSi while assuming a con- 
stant hopping integral t(i) — t and an elastic energy cost 
^KJ2 i S'f. This effect has been studied numerically in 
a different context jllj but it is relatively small in the 
case of the organic systems here considered. In contrast 
to the above on-site deformation, the positions of the 
intercalated anions can couple strongly to the electrons 
especially through modulations of the single particle hop- 
pings along the chains of the Peierls type, 



t(i)=t(l + Sf) 
with an elastic energy 

H ph = H elas = )-K B Y^(6; 
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The electron-lattice couplings have been absorbed in the 
re-definition of the displacements Sf so that the strength 
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of the lattice coupling scales like 1 /Kb (also t is set to 
1). We shall first consider the case of the isolated Hub- 
bard chain and, next, the role of an inter-chain magnetic 
coupling (in mean field). 

In contrast to previous treatments ]To| our numerical 
method enable us to obtain the lowest energy equilibrium 
lattice configuration without making any assumption on 
the broken symmetry ground state (GS). In particular, 
no super-cell order is imposed a priori and the GS config- 
uration is obtained through a self-consistent procedure. 
Indeed, the total energy functionals E({5f}) can be min- 
imized with respect to the sets of distortions {Sf} by 
solving the following non-linear coupled equations, 

K B 8f + t(cl a c l+Ua + h.c.) = . (4) 

Here (...) is the GS mean value obtained by ED (us- 
ing the Lanczos algorithm) of Hamiltonian (||) on cyclic 
L-site rings (with L up to 16 sites). Since the second 
term depends implicitly on the distortion pattern {8f}, 
Eq. (f|) can be solved by a regular iterative procedure 
A similar approach has been applied to the case 
of the on-site Holstein coupling O]. We should stress 
here that within the intrinsic limitations of the method 
(adiabatic lattice and finiteness of the system) our resolu- 
tion of the problem is basically exact (numerical accuracy 
better than 10 -7 ). Note also that once a small adiabatic 
lattice coupling is included finite size effects become quite 
small. 

Before discussing our main results on the Hubbard- 
Peierls chain, let us briefly describe the phase diagram 
of the Hubbard-Holstein chain |ll[] in order to introduce 
the generic types of CDW states. At quarter-filling, the 
Fermi wave vector is given by q2k F = f so that, at small 
U, one expects an instability towards a 2kp CDW state of 
wavevector X2k F — 4a (a is the lattice spacing) mediated 
by the electron-lattice coupling. In contrast, for large U, 
the system becomes more similar to a gas of interacting 
spinless fermions (SF) and the instability is likely to oc- 
cur at wavevector 2kp F — 4kp. More generally, we can 
parameterize the relative charge density modulation as, 

/\ t7 y ■ j° 

= P4k F COS (27T- 1 ) + p 2 k F COS (27T- L - + $2fc F ) 1 ( 5 ) 

n Za 4a 

where An^ = (riij — n. Complete phase diagrams of the 
Hubbard-Holstein chain have been established in Ref. p"l| ] 
and we briefly summarize them here. For V — the 
metallic uniform U phase (p2k F = Pik F = 0) is restricted 
to a region at small lattice coupling. Above a criti- 
cal line (1/K)u, three different insulating CDW phases 
can be distinguished; (i) at small U, a 2kp CDW phase 
{pAk F — 0) centered on the sites, i.e. with $2k F = 0; 
(ii) at intermediate U (in the range 4-8), a bond- centered 
2k F CDW phase (i.e. with $ 2 k F = f ); (hi) at large U, 
a Ak F CDW {p2k F — 0). As discussed in Ref. Q, a 
small NN repulsion suppresses completely the intermedi- 
ate phase and enlarges the region of stability of the Akp 
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FIG. 1. Typical (U/t,l/ Kb) phase diagrams of a ^-filled 
Hubbard-Peierls chain without (a) or with a nearest neighbor 
Coulomb repulsion V (b) obtained from ED of small periodic 
chains. Hashed regions are unphysical. 

CDW phase. Although these CDW might have some rel- 
evance to the low temperature phase of the (TMTTF^X 
family, their charge modulations should couple strongly 
to the anion potential. Therefore, we investigate next the 
role of a Peierls coupling. 

The phase diagrams as a function of the on-site re- 
pulsion U/t and the Peierls coupling 1 /Kb are shown 
in Figs. |l|(a,b) for V = and V — 2. For too large 
electron-lattice coupling the linear coupling approxima- 
tion breaks down and our model becomes unphysical 
(hashed regions) so that we shall restrict to small and 
intermediate values of 1/Kb- For intermediate values of 
I/Kb, the uniform state is unstable towards translation 
symmetry broken states. Such bond order wave (BOW) 
states are characterized by a modulation of the hopping 
amplitudes of the form, 

5? = 6f kF cos (2^£) + 6i kp cos (2tt^ + ^ kp ) . (6) 

Generically, we find that BOW coexist with weaker 
charge modulations (CDW) given by (|). The ampli- 
tudes of the bond and charge modulations are shown in 
Fig. ||(a) for a fixed electrondattice coupling Kb = 0.8. 
Very small finite size effects are observed and calcula- 
tions on 12- and 16-site rings give almost identical re- 
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suits. Two different types of structures are stable; (i) in 
the weak coupling regime (let's say U/t < 3), a strong 
2kp BOW with Q 2kp = i- e - corresponding to a X- 
X-Y-Y type of sequence of the bonds occurs. This mod- 
ulation coexists with a weaker 2kp site-centered CDW 
(A-B-A-B type of sequence of the on-site charge densi- 
ties, A and A corresponding to opposite values of Am) 
and an even weaker 4kp (i.e. A-B-A-B) CDW com- 
ponent (£>i phase); (ii) at larger U/t, the D 2 phase 
corresponds to the superposition of a lattice dimeriza- 
tion (4kp BOW) together with a tetramerization (2kp 
BOW with Q>2k =0)- I n other words, among the, let's 
say, weak bonds of the dimerized state, one every two 
becomes weaker (or stronger) so that electrons become 
weakly bound in singlet pairs on next NN bonds. The 
D 2 phase is therefore a simple realization of the above 
mentioned SP phase Jl3|. Interestingly enough, we ob- 
serve that the tetramerization leads to a weak coexisting 
2kp CDW component correspond to a A-A-B-B charge 
modulation (i.e. with &2k F — 7r /4). Note also that the 
boundary between the Z?i and D 2 phases is a first-order 
transition line as it is clear from the discontinuity of the 
various order parameters seen in Fig. ^(a). We have also 
included in this plot the energy difference between the 
states corresponding to DI and D2 patterns. This differ- 
ence is in general (except at large U) considerably smaller 
than charge and spin excitations (see below). 

To complete our study we have also computed the 
charge (A p ) and the spin (As) gaps in the D\ and 
D 2 phases for the same set of parameters as shown in 
Fig. ^(b). As for the Fourier amplitudes, finite size ef- 
fects are almost negligible especially for As- Clearly As 
follows closely the magnitude of the 2k F BOW-CDW. 
For large-U (and large dimerization) the system behaves 
qualitatively like a spin-1/2 antiferromagnet (since the 
electrons are localized on the strong bonds) and the spin 
gap is expected to vanish in this limit. In contrast, in 
the Di phase, electrons are strongly localized in pairs on 
2 adjacent strong bonds (i.e. on 3 sites) so that As ~ t. 
The charge gap, also shown in Fig. 0, has a minimum at 
intermediate U in the region corresponding to the cross- 
over from dominant 2kp to dominant 4kp BOW-CDW. 
Note that both charge and spin gaps are discontinuous 
at the first order transition between D\ and D 2 . 

Our study shows that the electronic correlations (both 
U and V) are essential to stabilize the D 2 phase which, 
we believe, is a fair realization of the SP phase of the 
(TMTTF) 2 X material. Note that the 2k F (4k F ) Fourier 
components are suppressed (increased) as U ft increases 
as seen in Fig. §(a) and in agreement with previous nu- 
merical calculations of the on-site charge-density and NN 
charge-transfer response functions of the extended Hub- 
bard model 0. 

Lastly, we investigate the role of an inter-chain cou- 
pling which is relevant e.g. in the case of (TMTTF^Br 
(or (TMTTF) 2 PF 6 under pressure). These systems 
which are less anisotropic than (TMTTF^PFg at am- 
bient pressure exhibit an antiferromagnetic (AF) phase 
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FIG. 2. (a) Amplitudes of the CDW (p2k F and Pik F ) 
and BOW (<5^, and Sf kp ) components versus U/t for 
(Kb)' 1 = 1.25 and V = 2 (computed on a 16-site ring). Open 
(close) symbols correspond to CDW (BOW). Note that, al- 
though identical symbols have been used, the 2k_F orders are 
in fact different, i.e. have different phases 3>2fc F and $|fc F i n 
the Di and D2 phases (see text). The crosses indicate the 
energy difference (in absolute value) between the states with 
these two patterns, (b) Charge (A p ) and spin (As) gaps in 
unit of t vs U/t computed on 12- and 16-site rings (As on 
L — 12 and L = 16 are undistinguishable). 

at low temperature. In an insulating regime, due to the 
presence of a charge gap A p , the inter-chain single parti- 
cle hopping t± is believed to be irrelevant ||. Therefore, 
we shall only consider a transverse magnetic coupling Jj_ 
(typically J± ~ tj_/A p ) in mean field approximation. 
Our previous method can be straightforwardly extended 
to include this inter-chain coupling by adding to ([!]) an 
extra term as, 

H x = i Yl H ^ n ^ - n *;l)' ( 7 ) 

i 

where the local fields Hi are determined self-consistently 
(simultaneously with the 8f bond modulations) from an 
additional set of non-linear equations, 

fli = ^(KT-niii)>- (8) 

Our results on coupled Hubbard-Peierls chains are 
summarized in Fig. Besides the D 2 phase (Fig. |i|(a)) 
which is stable at small Jj_ new magnetic phases de- 
picted in Figs. (|(b,c) appear; (i) at small electron-lattice 
coupling an antiferromagnetic phase consisting of a site- 
centered 4ki? CDW with a finite spin density on the sites 
carrying an excess charge (see Fig. f|(c)); (ii) at larger 
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FIG. 3. Phase diagram of weakly coupled i-filled Hub- 
bard-Peierls chains as a function of 1/Kb and Jj_ for [7=6, 
V — 2. Hashed regions are unphysical. 



coupling and the Coulomb repulsion. We have shown 
that, generically, lattice modulations (BOW) are always 
accompanied by CDW's of weaker amplitudes. In ad- 
dition, at intermediate and large on-site Coulomb repul- 
sion, the lattice modulation consists of a superposition of 
a 4fcp (dimerization) and a 2Uf (tetramerization) BOW. 
Interestingly enough, we found that a NN electronic re- 
pulsion further stabilizes this lattice/charge modulated 
phase. Under the application of an inter-chain AF cou- 
pling, we found new long range spin order phases showing 
coexisting 4kp CDW or BOW (i.e. dimerization). It is 
argued that such a simple model can well describe the 
various low temperature SP, AF and SDW phases of the 
insulating charge transfer salts of the sulfur series Jj| . 

Computations were performed at the Supercomputer 
Computations Research Institute (SCRI) and at the Aca- 
demic Computing and Network Services at Tallahassee 
(Florida) and at IDRIS, Orsay (France). Support from 
ECOS-SECyT A97E05 is also acknowledged. 
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FIG. 4. Pictures of the various phases of antiferromagneti- 
cally coupled chains showing the bond modulation (thick lines 
are strong bonds, dashed lines are weak bonds), the charge 
modulation (full and open circles correspond to excess and 
depression of charge respectively) and the local spin densities 
(arrows), (a) dimerized/tretramerized D2 phase; (b) Spin 
density wave SDWi phase; (c) antiferromagnetic phase. 



values of I/Kb, a superposition of a dimerization (4kp 
BOW) with a 2k F bond-centered SDW order |J| (see 
Fig. [|(b)). Note that only the spin densities of these two 
magnetic phases could be obtained by adding two out-of- 
phase 2kp CDW for the spins up and down so that the 
4kp CDW and BOW components should really be consid- 
ered as extra coexisting orders. In a small region of the 
parameter space, a more exotic magnetic SDW2 phase 
(not shown) has been stabilized on our clusters. This 
phase contains all three CDW, SDW and BOW com- 
ponents with large supercells of the order of our chain 
lengths. Finally, we note that the region of stability of 
the D2 phase should be extended by a small interchain 
bond coupling Jul. 

To summarize, the role of Peierls electron-lattice cou- 
plings has been investigated in the adiabatic approxi- 
mation in quarter-filled isolated or weakly coupled one- 
dimensional Hubbard chains. A numerical method based 
on ED techniques supplemented by a self-consistent pro- 
cedure has been used to determine the various phase 
diagrams as a function of the strengths of the lattice 
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